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Стохастическое имитационное моделирование механизмов укорочения 
теломер клеток в процессах старения и развития патологических отклонений” 


А. А. Бутов, М. А. Карев, С. А. Хрусталёв 


Рассматриваются вопросы применения стохастического имитационного моделирования к описанию и иссле- 
дованию механизмов укорочения концевых структур ДНК — теломер. Целью настоящей работы является 
исследование на основе стохастического и имитационного моделирования протекающих в клетках явлений, 
которые приводят к уменьшению длины теломер и, как следствие, к снижению пролиферативног о потенциа- 
ла и возникновению патолог: ИЙ, в том числе онкологических. Математическая модель описывается в семи- 
мартингальных терминах. Анализ результатов моделирования основан на сопоставлении со статистическими 
данными биологических экспериментов. Построенные распределения клеток по длинам их теломер сопо- 
ставлялись с экспериментальными данными, полученными на культурах фибробластов человека с ИСПОЛЬЗО- 
ванием метрики Леви—Прохорова. Результаты работы могут найти применение в медико-биологических ис- 
следованиях старения и возникновения различных патолог: ИЙ, а также при решении ряда вопросов в области 
геронтологии. 
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Введение. В последние годы выявлены неизученные ранее механизмы старения и процессов 
опухолевой трансформации клеток. Этому способствовало развитие современных биотехнологий, 
а также большой интерес к изучению роли теломер и фермента теломеразы в клеточном цикле. 
Раскрыты механизмы, связанные с теломерами, которые служат причиной канцерогенеза. Тело- 
меры — участки ДНК на концах хромосом, состоящие из одинаковых нуклеотидных повторов. 
Главная роль теломер в клетке связана с их участием в процессе клеточного деления. Деление 
клетки сопровождается потерей части концевых участков ДНК. При этом возникает так называе- 
мая проблема «недорепликации» концевых участков ДНК — теломер [1]. Но кроме участия в про- 
цессах деления клетки теломеры защищают геном, образуя защитный нуклеопротеидный ком- 
плекс на концах хромосом. Этот комплекс стабилизирует хромосомы и предохраняет их от разру- 
шительных воздействий — таких, например, как деградация и слипание концевых участков хро- 
мосом, повреждение ДНК на концах хромосом, нарушение рекомбинации [1, 2]. В случае повре- 
ждения защитного нуклеопротеидного комплекса, образованного теломерами, хромосомы начи- 
нают сливаться с большей частотой, что в результате приводит к тяжелым генетическим анома- 
лиям, к хромосомным перестройкам и деградации хромосом. Кроме того, изменения или наруше- 
ния в структуре теломер могут затрагивать не только их собственные функции, но и экспрессию 
жизненно важных генов, находящихся в других районах хромосом [3]. 

Восстановление теломерных участков хромосом осуществляет фермент теломераза. Бла- 
годаря этому механизму хромосомная ДНК и связанные с ней гены защищены от повреждения во 
время деления клетки. Теломераза — специфический клеточный рибонуклеопротеин, который по 
механизму реакции является обратной транскриптазой. Теломераза принимает участие в синтезе 
активных центров для связывания теломерной ДНК с теломерными защитными белками [4]. В эм- 
бриональных и стволовых клетках благодаря экспрессии теломеразы поддерживается постоянный 
уровень длины теломер. Кроме того, повышенные уровни теломеразы содержатся в большинстве 


* Работа выполнена при поддержке Министерства образования и науки РФ (НИР, проводимые в рамках государственного 
задания Министерства образования и науки РФ на 2013 г., Программа стратегического развития УлГУ на 2012-2016 гг., 
Программа развития деятельности студенческих объединений УЛГУ на 2012-2013 гг.). 
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опухолей человека. В соматических клетках уровень теломеразы постепенно снижается, и, как 
следствие, происходит уменьшение длины теломер. 

В настоящее время возникновение и развитие многих патологий в организме человека, в 
том числе онкологических, связывают с нарушениями механизмов регулирования длины теломер 
[5]. Поэтому исследования данных механизмов, выявление зависимостей между характеристиками 
развития заболеваний и динамикой изменений теломерных участков ДНК — актуальные и важные 
задачи. Эффективным инструментом описания и изучения таких задач является математическое и 
имитационное моделирование. 

Цель настоящей работы — исследование протекающих в клетках явлений, которые приво- 
дят к уменьшению длины теломер. Как следствие, снижается пролиферативный потенциал и воз- 
никают различные патологии, в том числе онкологические. На основе настоящего исследования 
разрабатывается математическая имитационная модель. Она позволяет не только вскрыть и каче- 
ственно объяснить механизмы развития данных патологий, но и рассчитать их количественно — 
на основе анализа стохастических деградационных процессов теломер. 

Важным отличием настоящей работы от указанных выше является применение семимар- 
тингальных методов и теории случайных процессов при построении математической модели ме- 
ханизмов укорочения теломер клеток. Это позволяет наиболее адекватно отражать стохастиче- 
скую природу явлений, протекающих в клетках. Компьютерные имитационные реализации данной 
модели дают возможность обнаружить различные зависимости между характеристиками процес- 
сов, протекающих в клетках, и патологическими явлениями в исследуемых организмах. 
Основная часть. Анализ последних исследований. Моделирование клеточных циклов и про- 
цессов, приводящих к старению клеток и потере ими способности к пролиферации, изучалось во 
многих работах [6-16]. Рассматривались различные подходы к отражению протекающих в клетках 
явлений с помощью математических методов и их компьютерной реализации. 

Комплексный подход к рассмотрению механизмов укорочения теломер клеток представлен 
в статьях [6, 7]. Авторы этих работ предполагают, что теломеры клеток укорачиваются не только 
при делении, но и вследствие оксидативного стресса. Механизм увеличения скорости укорочения 
теломер, возникающий при повреждении митохондрий клеток, рассмотрен в [7]. В [6] в качестве 
дополнительной причины уменьшения длины теломер представлено накопление однонитевых 
повреждений теломер свободными радикалами. 

Модель снижения пролиферативного потенциала культуры клеток при развитии опреде- 
ляется вероятностями изменения состояний клеток (как способных к делению, так и неспособ- 
ных). В результате моделирования для каждого поколения клеток в некоторых работах получены 
соотношения чисел делящихся и не способных к пролиферации клеток. Однако в данной модели 
не рассматривается природа явлений, которыми обусловлена потеря клетками способности к де- 
лению. Похожие результаты дают модели, рассмотренные в работах [8, 9]. Методы моделирова- 
ния основаны здесь на ветвящихся процессах, что не всегда адекватно отражает природу процес- 
сов, протекающих в клетках. 

Марковские ветвящиеся процессы также были использованы при построении математиче- 
ской модели укорочения теломер в [10]. Здесь получены зависимости соотношения числа деля- 
щихся клеток, а также средней длины теломер от достигаемого популяцией РО (РоршаНоп дочЫ!д). 

В работе [11] представлена модель, основанная на предположении, что основной причи- 
ной снижения уровня пролиферативной активности клеток является укорочение их теломер. При- 
чем полагается, что не все теломеры клетки влияют на её способность к делению, а только их 
часть. К сожалению, в этой модели предполагается единственный механизм укорочения теломер, 
возникающий при делении клеток. Однако в реальных культурах наблюдаются и другие причины 
изменения длины теломер [12]. 
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Математическая модель механизмов поддержания длины теломер, а следовательно, и 
пролиферативного потенциала клеток за счёт комбинации неравного обмена теломер у сестрин- 
ских хроматид рассмотрена в работе [13]. 

В [14] рассматривается динамика укорочения длин теломер в клетках, объясняющая ста- 
рение организмов, синдром Вернера и бессмертие клеток, обусловленное ферментом теломера- 
зой. При математическом моделировании здесь использованы численные методы Монте-Карло. 
Сравниваются математические модели четырёх возможных механизмов укорочения теломер. Про- 
стейшая модель предполагает, что фиксированный участок теломеры теряется при каждой ре- 
пликации. Вторая модель основана на гипотезе, согласно которой динамика укорочения теломер 
зависит от их длины. В третьей модели величина укорочения теломеры фиксирована при каждом 
делении, но вероятность деления клетки зависит от длины теломеры. В четвёртой и укорочение 
теломеры, и вероятность деления зависят от длины теломеры. 

В основе работы голландских учёных [15] лежит гипотеза о прямой пропорциональной за- 
висимости интенсивности укорочения теломер от их длины. Здесь математическая модель пред- 
ставлена в терминах нелинейных дифференциальных уравнений, описывающих процессы укоро- 
чения теломер, динамика которых, в свою очередь, связана с пуассоновскими потоками. 

В рамках исследования учёного из США [16] проделана огромная работа по изучению раз- 

вития рака. В результате построена математическая модель процессов канцерогенеза, основан- 
ная на системе уравнений в частных производных. Она учитывает многие аспекты данного явле- 
ния, в том числе и механизмы укорочения теломер. 
Математическая модель. Математическая модель развития клеточной популяции построена на 
основе траекторных семимартингальных методов. Суть этих методов — описание закономерно- 
стей протекающих в клетках процессов в терминах стохастического исчисления. Факторы, влия- 
ющие на состояние клеток, но не известные биологу или не учитываемые при исследовании раз- 
вития клеток, представлены в таком описании в форме стохастических компонент случайных про- 
цессов. Это даёт возможность формулировать в строгом математическом виде биологические за- 
кономерности, определяющие динамику развития клеток. Данные закономерности в форме стоха- 
стических уравнений легко алгоритмизуются. Таким образом возникает возможность компьютер- 
ного имитационного моделирования процессов развития клеточных популяций. 

В настоящей модели процесс изменения числа клеток в популяции практически представ- 
ляет собой многомерный семимартингал. Изменения состояния клеток определяются скачками 
точечных (считающих) процессов. Такое представление позволяет моделировать независимые 
реализации процесса развития популяции. Также это даёт возможность в любой момент времени 
строить эмпирические функции распределения клеток по длинам их теломер и выявлять зависи- 
мости показателей популяции от параметров модели. Например, можно устанавливать, в частно- 
сти, следующие зависимости: 

— достигаемая популяцией клеток величина удвоений численности (РО) — от уровня свободных 
радикалов или интенсивности деления клеток; 

— время прекращения пролиферации популяции — от соотношения параметров модели; 

— относительное число делящихся клеток, а соответственно, и пролиферативный потенциал 
культуры — от РВ. 

В модели рассматриваются следующие механизмы изменения длины теломер: 

— регулярное укорочение теломер, возникающее при делении клетки; 
— укорочение вследствие накопления однонитевых повреждений теломер, вызванных оксидатив- 
ным стрессом. 

Также учитывается деление клеток с активной теломеразой, которое происходит без уко- 
рочения длины теломер вследствие наращивания концевых участков ДНК перед репликацией. 
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Предполагается, что популяция клеток состоит из двух групп. 

Т. Клетки с активной теломеразой. Деление происходит с последующим восстановлением 
теломерных участков ферментом теломеразой. 

П. Клетки, в которых отсутствует теломераза. При каждом делении происходит регуляр- 
ное уменьшение длины теломеры. По мере укорочения теломерных участков снижается интен- 
сивность деления клеток, а при достижении теломерами критической длины клетки утрачивают 
способность к пролиферации. 

Все клетки подвержены повреждениям ДНК свободными радикалами, вследствие чего 
также происходит уменьшение длины теломер. 

Количество клеток в каждой группе в момент времени Ё определяется функциями С (Е) и 


Е (Е) соответственно. Общее число клеток в популяции представляет собой сумму этих величин. 


Каждая группа клеток — это совокупность подгрупп, которые составляют клетки с различной 
длиной теломеры. Подобное представление даёт возможность наблюдать распределение популя- 
ции клеток по длинам теломер в любой момент времени. 


Пусть считающий процесс В° = (В представляет собой число актов деления клеток 
Е #>0 


при активной теломеразе (группа Г). Процесс В° = (ве ) — есть число клеток группы Т, у кото- 


рых произошло повреждение теломеры свободными радикалами. Процессы В” = (В: ) 


#20 ' 


К" = (& ре представляют собой число актов деления и повреждения теломер клеток в группе П. 


Число апоптозов клеток в группах Г и П определяются точечными процессами 02° = (2 ). си 
> 


р" =(5/)_‚ соответственно. 
120 
В ходе развития клеток снижается активность теломеразы (например, при дифференци- 
ровке стволовых клеток). Затем, по прохождении определённого числа актов деления (лимит 


Хейфлика) утрачивается возможность пролиферации. В рассматриваемой модели число клеток, 
утративших активность теломеразы (т. е. перешедших из Т во П группу), определяется считаю- 


щим процессом $ = (5, ),‚ (см. приложение). 

Таким образом, соотношения процессов В°, К°, В, К, 05, ГУ, $ определяют значения 
функций С (Е) и Е(Ё) (см. приложение). 

Все точечные (считающие) процессы, используемые в модели, являются неубывающими, 
поэтому представляют собой субмартингалы. Например, субмартингал В° = (ВР р согласно тео- 


Е 
реме Дуба — Мейера, допускает точное разложение В° = В° + М°, где Вг = | Ь°С (5) — ком- 
0 


пенсатор, возрастающий процесс ограниченной вариации, М° — мартингал. 
При моделировании непрерывного точечного процесса В° используется его дискретная 


модель, которая основана на следующем факте. Если компенсатор процесса В° представим в 


Е 
виде интеграла некоторой функции С (Е): В® = | 5°С ($) 45 ‚ то вероятность скачка процесса В° 
[4 Е 

0 


за малый промежуток времени ДЁ равна: 
РВ -В8 = 1} = 5°б (Е) АЁ+о(ДЕ), при ДЕ->0. 
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Поэтому процессы, построенные по последовательности случайных величин 


1 
ЕВ ы, : 1, свер. °С, 

=, 3,...: = 1! где р ' 
= 0, свер. 1-Л#°б у 


сходятся к распределениям процесса Вс при ДЕ ->0. 
При моделировании использовался шаг дискретизации ДЁ = 0,001. Случайные величины 
&, /[=1,2,3,... строились с помощью генератора псевдослучайных чисел с периодом 2". 
Величина 65°С (#) — интенсивность скачков (деления клеток) процесса В° . Это означает, 
что в момент Ё за малое время Л произойдет приблизительно °С (5) . А актов деления клеток. 


Интенсивности скачков остальных точечных процессов в модели определяются параметрами г°, 
Г, г’, 45, а* и $ соответственно. То есть для процесса Ю° интенсивность скачков имеет 
вид: г°С (Е), для процесса ВР — БЕ (Е) ит. д. 


Подобная конструкция позволяет построить модель эволюции культуры клеток с учётом 
рассматриваемых факторов, влияющих на изменение их состояний. Параметры модели имеют 
определённый биологический смысл и могут быть оценены в ходе анализа экспериментальных 


данных. Например, величины 5 и зе имеют смысл средней продолжительности клеточного 
цикла в группах Т и П соответственно. Параметры г° и г’ представляют собой интенсивность 
повреждения свободными радикалами теломер клеток в группах 1 и П соответственно. Параметры 
а° и а’ — средние скорости апоптоза клеток в группах Г и И соответственно. Величина $ — ин- 
тенсивность утраты клетками теломеразной активности. 

Обсуждение результатов моделирования. Компьютерное моделирование развития клеточ- 
ной популяции позволило проследить динамику распределений клеток по длинам их теломер. При 
этом рассмотрены параметры, которые соответствуют предположениям о приоритете различных 
факторов, определяющих характер распределений. Таким фактором может быть, во-первых, регу- 
лярное укорочение теломер из-за недорепликации концов ДНК в процессе деления. Во-вторых, это 
случайные разрывы, происходящие вследствие нерепарации свободнорадикальных повреждений 
или, например, воздействия ионизирующей радиации. Так, в ходе моделирования была установлена 
линейная зависимость средней длины теломер клеток от концентрации свободных радикалов в тка- 
ни (см. выше параметры гб и г”). Данный результат обобщает выводы, полученные в [17]. 

Заметим, что теломеры клеток человека и других млекопитающих укорачиваются в разных 
тканях с различной интенсивностью. Например, в печени при циррозах скорость укорочения воз- 
растает при усилении пролиферативных процессов. При этом необходимо учитывать оба фактора 
укорочения теломер, что позволяет в процессе моделирования воспроизвести особенности ре- 
ально наблюдаемых распределений. 

В некоторых тканях основную роль в укорочении теломер может играть неспособность 
клеток репарировать свободнорадикальные повреждения. Поэтому основным фактором, опреде- 
ляющим динамику распределений клеток по длинам их теломер, является механизм случайных 
разрывов. Примером такой популяции клеток, адекватно моделируемой предложенным методом, 
являются диплоидные фибробласты человека [18]. 

Построенные распределения клеток по длинам их теломер сопоставлялись с эксперимен- 
тальными данными, полученными на культурах фибробластов человека [18] с использованием 
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метрики Леви — Прохорова. Данная вероятностная метрика используется в качестве критерия 
достоверности выбранных параметров с заданным значением ошибки =*: 


= а ЕЕ нЕ [>0: Р®® (х -=)-=< Е" (х) < Е®® (х +=) +=, Ух ЕД} р (1) 
где Е“? (х) — функция распределения популяции клеток по длинам теломер, построенная на 


основе экспериментальных данных; Е" (х) — функция распределения, полученная в ре- 


зультате компьютерного моделирования. 

Использование данной метрики позволяет сравнивать эмпирические и модельные функ- 
ции распределения и фиксировать те параметры, при которых отличие одной функции распреде- 
ления от другой не будет превышать заранее заданную ошибку. 

На рис. 1 представлено качественное соответствие характеров моделируемых и реаль- 
ных экспериментальных распределений. Это свидетельствует о возможности применения и раз- 
вития данного подхода к моделированию воздействий различных факторов на популяции кле- 
ток — в частности, на раковые клетки при действии лекарственных препаратов или ионизиру- 
ющего излучения. 
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Рис. 1. Экспериментальная и модельная кривые распределений диплоидных фибробластов человека по длинам теломер 
(тысяч пар оснований — ТПО) для популяций различных возрастов 


В табл. 1 представлены начальные данные и время моделирования. Выбор размерности 
параметров моделирования зависит от типа исследуемых клеток. 

Итоги имитационного моделирования позволили установить, что данная модель явилась 
косвенным подтверждением физиологического объяснения предела Хейфлика. Определено сред- 
нее время достижения данного предела для популяции фибробластов человека. Выявлены зави- 
симости распределения времени прекращения пролиферативной активности популяции клеток от 
их распределения по длинам теломер. 

На самом деле предел Хейфлика трактуется в литературе достаточно примитивно, т. е. 
как определённое число делений клеток, после которого пролиферация прекращается. Однако, 
как показала и настоящая работа, развитие клеточной популяции представляет собой сложный 
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процесс, природа которого стохастична. Если рассматривать вероятностную суть данного явле- 
ния, то всегда присутствует некоторое распределение времени прекращения пролиферации кле- 
ток, т. е. её старения. Поэтому предел Хейфлика не чётко обозначен, и могут возникать различ- 
ные эффекты продолжительности пролиферативной истории клеток в зависимости от параметров 
клеточной популяции и условий её развития (скорость деления клеток, величина укорочения те- 
ломеры при делении клетки, концентрация свободных радикалов в ткани и т. д.). Эти особенности 
можно наблюдать с помощью имитационного моделирования и предложенной математической 
модели. Так, в ходе моделирования были получены зависимости скорости роста популяции от 
концентрации свободных радикалов в ткани и от средней длины теломер клеток. 






































Таблица 1 
Параметры имитационного моделирования процессов развития клеточной популяции 
Параметр Группа, обозначение РО10 РО10 РО49 РО55 
Ггруппа 9, (0) К=1,2,..., М(Ё) 4000 | 4000 | 4000 | 4000 
Начальное количество клеток 
П группа #, (0) А =1,2,..., М(Е) 4000 | 4000 | 4000 | 4000 
Тгруппа р° 0,158 0,151 0,161 0,163 
Интенсивность деления 
П группа 5” 0,1485 0,153 0,163 0,173 
Тгруппа а° 0,18 0,18 0,18 0,18 
Интенсивность гибели 
П группа а” 0,17 0,17 0,172 | 0,172 
Интенсивность свободнорадикальных повре- Тгруппа гб 104 10 10“ 10“ 
ждений п группа Г 10 10 10? 105 
Интенсивность утраты клетками теломерной ы 11 11 11 11 
активности 
Время т 63 77 87 94 




















Заключение. В настоящей работе рассматривались модели, применимые к описанию и исследо- 
ванию механизмов укорочения концевых структур ДНК — теломер. В качестве основы для моде- 
лей использовались стохастические дифференциальные уравнения. Цель настоящей работы — 
исследование на основе стохастического имитационного моделирования протекающих в клетках 
явлений, которые приводят к уменьшению длины теломер, следовательно, к снижению пролифе- 
ративного потенциала и возникновению различных патологий, в том числе онкологических. На 
основе данного исследования разработана математическая имитационная модель, которая позво- 
ляет не только вскрыть и качественно объяснить механизмы развития этих патологий, но и рас- 
считать их количественно — на основе анализа стохастических деградационных процессов тело- 
мер. Адекватность данной модели подтверждает сопоставление с экспериментальными данными, 
полученными на культурах фибробластов человека с использованием метрики Леви — Прохорова. 
Результаты работы могут найти применение в исследованиях старения и возникновения различ- 
ных патологий, а также при решении ряда вопросов в области геронтологии. 

Приложение. Подробное математическое описание модели. Пусть № (Е) — общее число 


клеток в момент времени ЕВ популяции, которая включает в себя клетки двух групп: 
М(Е)=6 (Е) +2 (6), (2) 
где С (Е) — число клеток с активной теломеразой; Е (Е) — число клеток, в которых данный фер- 


мент не действует. 
Число клеток популяции, способных к восполнению теломерных участков при делении за 
счёт активности теломеразы, выражается равенством: 


6), (2). (3) 
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Здесь М(Ё) — максимальная (по всем клеткам популяции) длина теломеры в момент времени & 
Она выражена в условных единицах, которые представляют собой среднее число пар оснований 
(ПО) нуклеотидов, теряемых вследствие «недорепликации» при делении клетки. Выражение 
9, (Е) — число клеток с активной теломеразой (группа Г), длина теломеры которых равна (. 


Значение процесса д, (Ё), К=1, 2, 3,..., М (Е) в каждый момент времени Ё выражается 
следующим соотношением: 


9 (Е) = 9х (0) + Вх (Е) + Хе, ()-5 к (Е) 0% (Е), К=1, 


]=К+1 
М(Е 


9% (Е) = 9 (0) + В» (Е) - ве (6+ о Кук (Е) 5% (Е) - Вр (Е), К=2М(И-Ь (9 


9, (= 9, (0) +88 (#)- а (®)-2° (©), К=М(#. 


1=1 





Опишем подробнее каждое слагаемое в (4). Начальное количество клеток группы Т с теломерой 
длиной К равно 9, (0). Процесс Вг (Ё), где К =1,2,..., М (Е), представляет собой число делений 
клеток с теломерой длиной К. Процесс Кр, (Е), где К =1,2,..., М(Ё), /=12,..., К -1, определя- 
ет число клеток, которые вследствие случайных повреждений свободными радикалами утратили 


участок теломеры длиной /. Поэтому общее число клеток, имевших теломеру ВИНО К, в кОТО- 


рых укорочение произошло из-за свободнорадикальных повреждений, равно Ус, ; (Е Соответ- 
1= 


ственно, число клеток, теломера которых имела длину, большую К, и уменьшилась до размера К 
м8 

вследствие случайных повреждений, равно »` К, , (Ё). 
1=К+1 


Процесс $, (Ё), К =1,2,..., М (Ё) представляет собой число клеток с активной теломера- 
зой и длиной теломеры К, которые в ходе развития утратили действие данного фермента (пере- 
шли в число клеток группы ТТ). Число апоптозов клеток группы Г определяется процессом 08 (Е) , 


Ким (Е 


Аналогично определяется число клеток, утративших теломеразную активность (группа 11): 
М(Е) 
Е -У,^ (1), (5) 
К=1 
где Ё, (Е) — число клеток, длина теломеры которых равна К. 
Соотношение, определяющее значения процесса #, (Ё) ‚ имеет вид: 


Е (Е) =Е (0) +28, (Е) + ;У К’ к (Е) +5, (#)- В (Е), К=Ь 


= 
К 


Г, (Е) =^ (0) В» (Е) +28... (Е) >.Вк, (В+ _ ее ВЕ а О) 


(Е) -Е (0)- В* () $40.52, к-М(е). 


Также приведём описание слагаемых в (6). Процесс В, (Е), где К =1,2,..., М (Е) , представляет 








собой число клеток, при делении которых происходит регулярное укорочение теломер двух обра- 
зовавшихся дочерних клеток на одну условную единицу. То есть длина теломеры этих клеток уже 
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равна К -1, и следовательно, выражение вычитается. Соответственно, число дочерних клеток, 
длина теломеры которых стала равной К после деления клеток с теломерой длиной К +1, равно 


Е 
28... (Ё). 
Предполагается, что все клетки подвержены свободнорадикальным повреждениям, по- 
этому смысл значений процессов К; ‚ аналогичен описанным выше процессам КР и 
Процесс $, (Е) также участвует в соотношении (4) и описан выше. Число апоптозов кле- 
ток группы П определяется процессом СД» (Ё), К =1,2,...,М(Ё). 


Распределения всех описанных выше стохастических процессов, участвующих в построе- 
нии математической модели, характеризуются своими компенсаторами”. Ниже приведены уравне- 
ния компенсаторов рассматриваемых процессов. 

Компенсатор процесса, определяющего число делений клеток группы Т, имеет вид: 


Е 
5с [2 
ВВ (Е) =5° [ 9, (Е). (7) 
0 
Здесь 5° >0 — интенсивность деления клеток, задаваемая как регулируемый параметр модели. 
1 ь ь 
Величина 55 имеет смысл среднеи продолжительности клеточного цикла клеток с активнои те- 


ломеразой. 
Компенсатор процесса, представляющего число случайных повреждений теломеры клеток, 
имеет вид: 


Кр, (Е) =1е, [ 9, (Е) а, (8) 


где Г. ‚ КЕ2,3,.... М (Е), Г=1,2,....К -1 — интенсивность свободнорадикальных повреждений 


теломер, которая пропорциональна концентрации свободных радикалов в клетках. 
Компенсатор процесса, представляющего число клеток, утративших теломеразную актив- 
ность, имеет вид: 


$, (Е) =$[ 9, (Ка, (9) 


где $ — интенсивность перехода клеток из группы 1 в группу П. 
Компенсатор процесса, считающего апоптозы клеток группы Т, равен: 


65 (в) = 9° [ 9, (#)а, (10) 


где 4° — интенсивность гибели клеток. 
Аналогичный вид имеют компенсаторы точечных процессов в соотношении (6), которое 


описывает динамику изменений клеток группы П. Соответствующие интенсивности (Б”, . ут а") 


также являются изменяемыми параметрами модели. 


* Для любого точечного (считающего) процесса, траектории которого представляют собой кусочно-постоянную функцию со 


скачками, равными 1, имеет место разложение: О, = О, +т,. Здесь т = (т,), — локально квадратично-интегрируемый 


Е 
мартингал, О = (©. й ‚ — компенсатор процесса О. Если компенсатор представим в виде О, = [9.45 ‚ то функционал 
в 0 


9 = (9,),„ представляет собой интенсивность скачков процесса @ , т.е. Р{0,,,-@, =1| } =4.-А+о(А), А->0. 
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ЗТОСНА$ТТС $ТМУГАТТОМ МОБЕММС ОР СЕН. ТЕГОМЕВЕ ЗНОВТЕМТМ 6 
МЕСНАМТ$М$ 1М АСЕТМС АМО ОТ$ТОВВАМСЕ ОЕУЕГОРМЕМТ РКОСЕ$$Е$” 


А. А. Вико, М. А. Кагем, $. А. КАги$а!еу 


5оте 155ие5 о! арр/утд {те юсвазНс зтшаНноп тоде/та 0 1е аезспрёоп апа зшаута оЕОМАЗ епа тис игез 
(1е/отегез) зпомепттд тесваттсх аге соп$!АЕгеа. ТПе айт о! 5 гезеагсй [5 Ю зшау — оп ве Базе оЁ 1е $ю- 
спазНс апа тиайоп тоде!тд — ргосеззез (п {Те се/5 [еа@ тд ю пе {еотеге5' епоён зропеп/пд, апа аз а соп- 
5едиепсе, {о десгеазтд {те ргоШегаёуе сарасйу, апа деуеортд ра поосбиез тс!иата опсообис опез. ТВе таЁй- 
етайса!/ тодЕ! 15 аезстрей т е зетита!тдее егтз. Тре апауя5 оЁ те тоде/та гези!з 15 Базе оп {Тег 
сотрайоп Ю {те Боосдгса! ехрейтепЕ5 аа. Тре Гогтеа се! @тривоп т ассогаапсе 0 те! еотегеэ” [епоЁТ 
5 сотрагеа ий! {те ехрептета! даа обатеа Епгоидн {те (е5Ё5 регогтеа ийт вВитап НгоБазЕ-си!иге итд 
[еуу — Ргоквогоу тес. Тпе $шоау гезий5 сап Бе изеа т Ыотесгса! гезеагсй оп адета, апа оп аеуеор/пд уап- 
диз райообиез, аз ие! аз п отд а питрег оЁ ргоМетз т {те Пе о! дегото/ду. 

Кеуигога5: чосразЁс тиайоп тоде/тд, зетйтатда!, сотрепзаёог, ОМА, (еотеге, {=отегазе, адта, саг- 
сподепе$6, охаёвуе °е5°, ее габса(. 


* Тне гезеагсН 15 зирромеа Бу е Митеу оЁ ЕдисаНоп апа 5аепсе о! {Не Виз$ап РедегаНоп (В&р сопаисе4 мт пе пате 
ОГ {Ле доуегптепЕ {азК оЁ ВЕ Мит&гу ог ЕдисаНоп апа 5аепсе Гог 2013, $ 5гаедс аеуеюртеп ргодгат Гог 2012—2016, 
Реуеюортеп* ргодгат оЁ ($ экидепЕ ипюп$ Рог 2012-2013). 


109 


